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AWARDS ABSTRACT 

. m»ii i-i scheme for estimsting the 

The invention relates to am g ^ carrie r signal phase 

parameters of a received carrier sig " * ’ bigh Doppler, Doppler rate, 

modulated by unknown data expenenc g ^ 8 Global positioning System 

etc., that arises, for example, in dlrectly related to the position, 

(GPS) where the signal parameters receiver relative to three or 

velocity, acceleration and jerk of the GPS receiver 

four satellites in the system. 

FIG. 1 illustrates a block diagram of a 

parameters estimation by diff er ® nt1 ^ ^^gigna^parameters estimation by 
FIG. 2(a) illustrates a block and FIG. 2(b) illus- 

least square (LS) signal 6 parameters estimation by least 

trates a conceptual block g noise . FIG. 2(c) illustrates a 

square (LS) algorithm with prew l ma ti 0n by least square in an ap- 

block diagram for signal parameters y pic . 3 ( a ) illustrates a 

proximate realizable equivalent of FI * # * fi t stage of the present 

generalized pirates a Actional block model for a more 

invention, and FIG. 3(b) eneralize d single-stage estimator shown in 

compact representation of the g . , bloC k model for subsequent 

FIG. 3(a). FIG. 3(0 Illustrated a functional blook^^ * ln 

stages of a multistage estima or 3 used in simulated analysis of 

graphs a, b and o, a m d J t X.e eawSr shown in FIG. 1 where „-2. 

tracking performance of a multistage estimauoi 

Respective FIGS. 6 through 15 are graphs of, 
losing frequency lock vs CNR for ' * erroI . 8 vs CNR for a DUS algorithm in 
modulation; RMS frequency es nrobability of losing frequency lock vs 

the absence of data sence 0 f data modulation; RMS frequency 

CNR for a DLS algorithm in the presence presence of data modu- 

estimation error vs. CNR for vs CNR for a DLS-EKF 

1 at i on ; the probability of losi g 9 frequency estimation error 

algorithm in the absence of data , no ^diation J.MS ™“^ lation . modulo 2 > 

vs CNR for a DLS-EKF algorithm in K f al eorlth m in the absence of data 
RMS phase error vs CNR for a vs CNR for a DLS-EKF 

modulation; the prooabili y modulation' RMS frequency estimation 

algorithm in the presence of data f " 0 ^ l3 “ 0 "’ eaenoe of data modulation; 

error vs CNR for a DLS ‘™ , a 8 ° ror w3 C NR for a DLS-EKF algorithm (with 
and RMS pseudo range estimation err 

and without data modulation). 

k of the more general multistage scheme, 

In a two-stage embodiment . t scma res algorithm shown in 

the first stage (selected to be a _ mo ■ l re3ulting q in hig her rms estimation 
FIG. 3a) operates as a 9 oa T® e Drobab ility of the frequency estimation 

errors but with a relatively amnl \ n(r frequency. It provides relatively 

coarse "estimates °of Se^requency and ^ its der ^ operates^ of the 

- overall esbimabes 




MULTISTAGE ESTIMATION OF RECEIVED 
CARRIER SIGNAL PARAMETERS UNDER VERY 
Hlffl dynamic CONDITIONS OF THE RECEIVER 


ori gin of the Invention 

The invention described herein was made in the perform- 
ance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-517 (35 USC 202) in which the 
5 Contractor has elected not to retain title. 


T echnical Field 

The invention relates to a multistage scheme for es 1 - 
mating the parameters of a received carrier signal, such as a 
10 carrier signal phase modulated hy unknown data experiencing 
very high Doppler, Doppler rate, etc., that arises, for exam- 
ple, in the case of a Global Positioning System (GPS) whei e 
the’ signal parameters are directly related to the position, 
velocity, acceleration and jerk of the GPS receiver relative 
15 to three or four satellites in the system. 

Back ground Art 

The problem of estimating the parameters of a receive 
quasi- sinusoidal signal in the presence of noise occurs In 
diverse scientific and engineering disciplines. The signal 
parameters of interest are usually the phase, frequency and 
frequency derivatives which are varying with time. The esti- 
mation problem becomes considerably more difficult if the 
received carrier is modulated by unknown data while slmultane- 
25 ously experiencing very high Doppler and Doppler rate. 

just noted, situations occur In the cases of Global Position- 
ing System (GPS) receivers, but may also occur in NASA deep 
space communication links under high spacecraft dynamics. 

in a paper by W.J. Hurd, J.I. Statman and V.A. Vlln- 
30 rotter titled "High Dynamic GPS Receiver Using Maximum Liken- 
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of the phase along with a more refined estimate of frequency as well and in 
the process also reduces the number of cycle slips. 

The novelty of the invention resides in using a first stage for 
coarse estimation in cascade with a second stage which is then subject to 
significantly less dynamic changes in its signal input, thereby accommodat- 
ing higher dynamics in the signal received due to high dynamics in the 
motion of the receiver platform. 
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hood Estimation and Frequency Tracking," IEEE Trans., Vol. 
AES-23, no. H, pp. 925-837, July 1987. an estimator scheme is 
proposed and analysed for the OPS system based on the maximum 
iikelihood estimation (MLE) of code delay and Doppler fre- 
quency over a single data bit period. This scheme estimates 
Doppler frequency (assumed constant over successive intervals 
of bit periods) and then determines frequency and frequency 
rate by a Kalman filter tracking Doppler. The scheme does not 
involve carrier phase estimation. For the dynamic trajec- 
tories simulated in the paper by W.J. Hurd, et al . , the 
approximate MLE performance exhibited a threshold of about 
dB-Hs in terms of the received carrier power-to-noise power 
spectral density ratio (CNR), below which rapid performance 
deterioration occurred. 

For CPS applications, an alternative scheme has been 
proposed by C.E. Hoefener and L. Wells in a paper titled 
"Utilizing GPS for ultra-High Dynamic Vehicle Tracking in 
Space " Proceedings of the International Telemetering Confer- 
ence, Las Vegas, pp. 771-773, October 1986. wherein a parallel 
(nondynamic) link is established between the GPS satellites 
and a control ground receiver for the purpose of communicating 
the data to the ground receiver. The ground receiver simulta 
neously receives the frequency translated version o! the GPS 
receiver signal and removes the data modulation from this 
dynamic signal. Such an effectively demodulated signal is then 
processed by the estimation algorithm to obtain the desired 
signal parameter estimates. There are several estimation 

schemes in the literature for this problem. See, for example, 
R. Kumar, "Fast Frequency Acquisition via Adaptive Least 
Squares Algorithm," Proceedings of the International Tele- 
metering Conference. Las Vegas, pp. 9.-101, October ,9861 R. 
Kumar, "Fast Frequency Acquisition via Adaptive Least Square 
Algorithm," IEE Proceedings, Vol. 136, Pt. F, No, 1. PP- 
,55-160, August 1989, R. Kumar, "Differential Sampling for 
i Fast Frequency Acquisition via Adaptive Extended Least Squares 
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Algorithm," Proceedings of the International Telemetering 
Conference, San Diego, pp. 191-201, October 1987. 

More recently in R. Kumar, "Efficient Detection and 
Signal Parameter Estimation with Applications to High Dynamic 
5 GPS Receiver," JPL Publication 88- 1 12, National Aeronautics and 
Space Administration, Jet Propulsion Laboratory, California 
Institute of Technology, Pasadena, California, a scheme for 
simultaneous detection and estimation has been proposed. This 
scheme is based upon first estimating the received signal's 
10 local (data dependent) parameters over two consecutive bit 
periods, followed by the detection of a possible jump in these 
parameters. The presence of the detected jump signifies a 
data transition which is then removed from the received sig- 
nal. This effectively demodulated signal is then processed to 
15 provide the estimates of the global (data independent) parame- 
ters of the signal related to the position, velocity, etc., of 
the receiver . From the simulations, it is seen that the 
scheme offers very significant improvement in terms of the 
required CNR over the approximate MLE algorithm of W.J. Hurd, 
20 et al . , (1 987) cited above. A key feature of this scheme, 
which has a computational complexity of about three times that 
of a single extended Kalman filter, is that to a certain 
extent the data detection is independent of the acquisition of 
the phase or frequency of the received carrier signal in 
25 contrast to the conventional decision-directed phase-locked 
loop receivers in which the data detector is an integral part 
of the loop and depends upon the acquisition of the carrier 
phase and/or frequency. Thus, under low CNR and/or high 
dynamic conditions, the loop may not acquire lock or frequent- 
30 ly lose it during tracking. 

An object of this invention is to provide an efficient 
method for estimating the parameters of a received carrier 
signal without undue computational complexity utilizing a 
multistage scheme. 
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Statement of the Invention 

In accordance with the present invention, the first 
stage of a multistage estimator operates as a coarse estimator 
resulting in higher rms estimation errors but with a rela- 
5 tively small probability of the frequency estimation error 
exceeding one half of the sampling frequency (an event termed 
cycle slip). The second stage of the multistage estimator 
operates on the error signal available from the first stage, 
refining the overall estimates, and in the process reducing 

10 the number of cycle slips. 

The first stage algorithm is preferably selected to be 
a modified least squares algorithm referred to as differential 
least squares (DLS) algorithm. This estimation stage provides 
relatively coarse estimates of the received signal frequency 
15 and its derivatives. The second stage algorithm is preferably 
a third-order extended Kalman filter (EKF) which yields a more 
refined estimate of frequency as well as an estimate of the 
signal phase. 

A major advantage of the proposed multistage estimator 
20 is a reduction in the threshold on received carrier power-to- 
noise power spectral density ration (CNR) as compared to the 
threshold achievable by either of the two cascaded estimator 
stages alone. In fact, it appears from simulations that for 
the case of an unmodulated carrier, the proposed scheme 
25 achieves the same threshold as for an almost exact and compu- 
tationally intensive implementation of the maximum likelihood 
estimator (MLE) . For the case of a data modulated carrier, 
the proposed scheme provides an improvement of about 6 dB in 
terms of CNR compared to an earlier approximate MLE scheme 
30 reported by W.J. Hurd, et al . , (1987) cited above. The over- 
all complexity of the algorithm is about two times the com- 
plexity of a third-order Kalman filter or a single 

fourth~order extended Kalman filter • 

The novel features that are considered characteristic 
35 of this invention are set forth with particularity in the 
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appended claims. The invention will best be understood from 
the following description when read in connection with the 
accompanying drawings. 

Brief Description of the Drawings 

FIG. 1 illustrates a block diagram of a prior-art 
system for signal parameters estimation by differential least 
squares (DLS) algorithm. 

FIG. 2(a) illustrates a block diagram for signal pa- 

rameters estimation by least square (LS) algorithm with model 
noise colored. 

FIG. 2(b) illustrates a conceptual block diagram for 

signal parameters estimation by least square (LS) algorithm 

with prewhitened noise . 

FIG. 2(c) illustrates a block diagram for signal pa- 

rameters estimation by least square in an approximate realiz- 
able equivalent of FIG. 2(b). 

FIG. 3(a) illustrates a generalized single-stage estim- 
ator used in the first stage of the present invention. 

FIG. 3(b) illustrates a functional block model for a 

more compact representation of the generalized single-stage 
estimator shown in FIG. 3(a). 

FIG. 3(c) illustrates a functional block model for 

subsequent stages of a multistage estimator. 

FIG. *1 illustrates a multistage estimator for the 
process Gj(t). 

FIG. 5 illustrates in graphs a, b and c, a high dynamic 
trajectory used in simulated analysis of tracking performance 
of a multistage estimator shown in FIG. 4 where n=2. 

FIG. 6 is a graph of the probability of losing fre- 
quency lock vs CNR for a DLS algorithm in the absence of data 

modulation. 

FIG. 7 is a graph of RMS frequency estimation error vs 
CNR for a DLS algorithm in the absence of data modulation. 
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FIG. 8 is a graph of the probability of losing fre- 
quency lock vs CNR for a DLS algorithm in the presence of data 
modulation . 

FIG. 9 is a graph of RMS frequency estimation error vs. 
CNR for a DLS algorithm in the presence of data modulation. 

FIG. 10 is a graph of the probability of losing fre- 
quency lock vs CNR for a DLS-EKF algorithm in the absence of 
data modulation. 

FIG. 11 is a graph of RMS frequency estimation error vs 
CNR for a DLS-EKF algorithm in the absence of data modulation. 

FIG. 12 is a graph of modulo 2it RMS phase error vs CNR 
for a DLS-EKF algorithm in the absence of data modulation. 

FIG. 13 is a graph of the probability of losing fre- 
quency lock vs CNR for a DLS-EKF algorithm in the presence of 
data modulation. 

FIG. 14 is a graph of RMS frequency estimation error vs 
CNR for a DLS-EKF algorithm in the presence of data modula- 
tion. 

FIG. 15 is a graph of RMS pseudo range estimation error 
vs CNR for a DLS-EKF algorithm (with and without data modula- 
tion). 

Detailed Description of the Invention 

1 . Introduction 

An alternative scheme for the estimation of the signal 
parameters will now be described for both the case of unmodu- 
lated carrier signal and the case of a carrier signal phase 
modulated by unknown data. The proposed scheme is somewhat 
simpler than that of R. Kumar, described in the aforesaid JPL 
Publication 88-42, in that it is not essential to explicitly 
detect the data modulation for the second case. Basically, 
the new algorithm involves an appropriate modification of the 
DLS scheme of R. Kumar in the aforesaid paper titled "Differ- 
ential Sampling for Fast Frequency Acquisition via Adaptive 
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Extended Least Squares Algorithm" so as to apply the algorithm 
to the case of unknown data modulation. 

As discussed in that paper, if a DLS technique is 
applied with the Nyquist sampling of the received signal, a 
loss in performance is expected compared to the optimum 
achievable performance. Consequently, oversampling and cyclic 
sampling was proposed to avoid such a loss. In the present 
invention, sampling at Nyquist rate is proposed with an alter- 
native method for estimating parameters in order to keep the 
overall performance close to optimum. The new scheme proposed 
consists of a multistage procedure wherein the parameters of 
the signal are estimated in more than one stage. First, coarse 
parameters are estimated by an algorithm like DLS which has a 
low threshold on CNR but with possibly higher rms estimation 
errors. Then an error signal whose parameters are equal to 
the difference between the true parameters and the coarse 
estimates is processed by another algorithm to estimate these 
error signal parameters. Since the error signal involves much 
smaller dynamics, the second algorithm can have smaller band- 
width resulting in a smaller estimation error. In principle, 
this procedure of processing an error signal in another stage 
can be repeated any number of times with successive stages 
having progressi vely lower bandwidths. 

The example described below of a multistage estimator 
is confined to two stages of recursion and applies a modified 
least square algorithm for the first stage and a third-order 
extended Kalman filter algorithm for the second stage. It was 
expected that the overall algorithm would have both smaller 
threshold and smaller estimation errors compared to either 
algorithm operating by itself. Indeed, this is borne out by 
simulations presented hereinafter. Thus, for the case of no 
data modulation, the threshold on SNR is about 1.5 dB lower 
than the third-order EKF , and the estimation errors are only 
marginally higher than for the third-order EKF alone. The 
threshold achieved is in fact the same as achieved for a 
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nearly exact implementation of the maximum likelihood estim- 
ator (MLE). It is also noted that the threshold achieved by 
the present invention is about 3 dB lower than conventional 
cross product AFC (CPAFC) loops and phase-locked loops com- 
5 pared by V. A. Vilnrotter, S. Hinedi and R. Kumar, "A Cornpari 
son of Frequency Estimation Techniques for High Dynamic 
Trajectories," JPL 88-19, Jet Propulsion Laboratory, Califor- 
nia Institute of Technology, Pasadena, California, whereas the 
rms error is less than one half of that obtained by CPAFC. In 
10 the simulation described hereinafter , the rms error is margin 
ally higher than for a third-order EKF due to the nonoptimal 
sampling used in the DLS algorithm. 

For the case of data modulation, results of the simula- 
tions are compared with those reported by Hurd, et al . , 
15 (1987), where analysis and simulations are presented on the 

performance of Fast Fourier Transform ( FFT ) based MLE algo- 
rithm. In that report, the trajectories of the GPS signals 
have somewhat less severe dynamics compared to those consid- 
ered in this presentation. In terms of threshold on CNR, the 
20 proposed scheme of this presentation exhibits a threshold of 
24 db-Hz compared to about 30 dB-Hz reported by Hurd, et al . , 
thus providing an improvement of about 6 dB. In terms of the 
rms frequency estimation errors at a 30 db-Hz CNR, the scheme 
of Hurd, et al . , provides a rms range rate error of about 6 
25 m/ s compared to an error of less than 2 m/s achieved in this 
simulation. There is also a very significant improvement in 
terms of the rms position estimation error. At about 30 
dB-Hz, an rms error of 1 meter is reported by Hurd, et al . , 
compared to about 0.25 meter obtained by the proposed algo- 
30 rithm. It may also be remarked that in the previous scheme of 
Hurd, et al . , pseudo-random codes with rate 10.23 MHz are 
needed for the purpose of range measurements, thus requiring a 
zero-crossing channel bandwidth in excess of 20 MHz. The 
present invention on the other hand, extracts the range inf or 
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mation from the carrier signal itself and thus needs a band- 
width equal to only a fraction of 1 MHz. 

2. Receiver Configuration 

Consider the problem of estimating the high dynamic 
phase process OjCt) of the desired signal Sj(t) observed in 
the presence of an additive narrow-band noise process v j- ( t ) as 

TjCt) = SjUJ+VjCt) 

= ASin(uj c t+0 I (t)+TrD(t))+v I (t), (1) 

where u> c is the received signal carrier frequency in the 
absence of any dynamics and D(t) is a binary digital waveform. 
In the case of a GPS receiver, the dynamic phase process 0j(t) 
arises from the receiver dynamics and, over a sufficiently 
small estimation period, 

0 : (t) = 0 Io ^ Io t 1 + Y Io t 2 + 6 Io t 3 (2) 

for some unknown parameter vector * lQ = [e Io w l0 Y Jo 6 l0 ]. In 
a somewhat simpler version of the problem, the data modulation 
D(t) is either absent or is assumed known and thus can be 
eliminated from Equation (1). Both of these cases will be 
treated in some detail. 

In the first stage, the present invention estimates the 
parameters related to the frequency and its derivatives using 
the DLS algorithm. For this purpose, the received signal 
Sj(t) is quadrature demodulated in a conventional input sec- 
tion 10 comprised of mixers 11 and 12 , a ir/2 phase delay 
element 13 which shifts the phase of a signal s L (t) from a 
voltage-controlled oscillator (VCO) 14 and lowpass filters 15 
and 16 as shown in FIG. 1. The VCO frequency update rate of 
T/Tg indicated by a block 17 is controlled by a signal 0(N) 
produced by a differential least square (DLS) algorithm indi- 
cated by a block 18 the function of which will be described 
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below. Assuming that the input to the VCO is a signal 0 L (t) 
which is an appropriate quadratic function of time, t, result- 
ing in the following VCO output signal s L (t) 

5 s L (t) = 2Cos(u) c t+0 L (t)) 

e L (t) - 0 Lo +w Lo t+T Lo t 2 + 6 Lo t3 (3) 

for some constant vector 4 i Lo = [0 Lo <*> Lo T Lo S Lo ]', the sam- 
10 pled version of the in-phase and quadrature components of the 
demodulated signal are given by 

y (k ) = ASin(e(k) + irD(k)) + v^k) 

15 z(k) = ACos(0(k) + nD(k)) + v q (k) ;k = 1,2 N (4) 

where 


0 ( k ) = QjOO - 0 L (k) = 0o^o kT s +Y o (kT s )2+6 o (kT s )3 

20 

4*0 = ♦io “ ^Lo = ^ 3 * * 6 o w o T o 6 o ] ' 

and ij> 0 i 3 the parameter vector characterizing the error signal 
to be estimated, with T s denoting the sampling interval. In 
25 Equation (4), v^k) and v q (k) represent the sampled in-phase 
and quadrature components of the bandpass noise process Vj(k). 
The parameter vector is estimated by the differential least 
squares ( DLS) algorithm of R. Kumar (1988) cited above, as 
described in the following section. 

30 

3. Differential Least Squares (DLS) Algorithm 

Consider first the problem of estimating the unknown 

parameters u 0 , Y 0 and 6 0 from the measurement Equation (4) for 

the case of no data modulation, i.e., when D(k) = 0, and 
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expand Sin(0(t)) in a Taylor series around t k _-j = (k-1)T s to 
obtain 


Sin(G( t ) ) = Sin (0(k-1 ) ) +T0(k-1 )Cos(0(k-1 ) ) +. . . (5) 

5 with a similar expansion for Cos (Q(t)). For small (t-t (< _ 1 ), 
the series in Equation (5) may be approximated by the first 
two terms and from Equation (iJ) the following differential 
signal model with t k = (k-(1/2))T : 

10 y d (k) £ y(k)-y(k-1)=T s (a) 0 +2Y 0 T k +36 0 T^)z(k-1)+C i (k) 

z d (k) ^ z(k)-z(k-1)=-T s (cj 0 +2Y 0 T k +3<5 0 T k )y(k-1)+c q (k) (6) 

where 

15 

q(k) J v i (k)-v i (k-1)-T s (m 0 +2Y 0 T k +36 0 T^)v q (k-1) 

?i(k) i v Q (k)-v q (k-1)-T s (u) 0 +2Y 0 T k +36 0 T^)v i (k-1) (7) 

20 The measurement model of Equation (6) may be rewritten in the 
following standard: 

Z d (k) » H' (k)B+£(k) ; k=1 , 2, • • . N (8) 

25 where 


30 


B' - U c 2 Y q 66 q ] 

’t z(k-1) T 3 T k z(k-1) 0.5T s i^z(k-1) " 

k > = 2 

_-T 3 y(k-1) -T 3 i k y(k-1) -0. 5T s T^y (k- 1 ) 


(9) 


z ^lO<) = [y d (k) z d (k)] ;C(k) = (k ) 5 q (k)] 
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The parameter vector 3 in Equation (8) is now estimated by an 
exponential data-weighted least squares algorithm in a recur- 
sive or nonrecursive form (Kalman filter). in its nonrecur- 
sive form, the estimate of 3 obtained on the basis of N 
5 measurements and denoted by 3(N) is given by 

N . N 

e(N) = { LH(J)H'(J)A M “J}’ 1 ( J H(j)Z d (j)A N "J} (10) 

J 1 J ® I 

where A is some appropriate weighting coefficient within 
10 0<A <1 . An equivalent recursive form of Equation (10) is the 

following algorithm: 

3(N) = 3(k-1 )+L(k)e (k) 

15 L(k) - P(k-1)H(k)[AI + H' (k)P(k-DH(k)] -1 

( 11 ) 

P(k) = (P(k-1)-P(k-1)H(k)[AI+H , (k)P(k-1)H(k)] _1 H'(k)P(k-1)}/A 
c(k) = Z d (k)-H'(k)3(k-1) ; k = 1,2 N 

20 

Note that the matrix to be inverted in Equation (11) is only a 
(2x2) matrix. In an alternative but equivalent form, one may 
process the scalar measurements y d (k), z d (k) sequentially 
instead of working with the vector measurement Z d (k). Moreo- 
25 ver , the matrix P(k) of Equation (11) with k=N is the same as 
the matrix inverse in Equation (10), i.e., 


P _1 (k) = .1 H(j)H'(j)A k "J 
j ' 


30 


Alternatively, the matrix P -1 (k) may be written as 


( 12 ) 


P'1 


(k) = l A 


35 


L A {z 2 ( j )+y 2 ( j ) } Bi T^ : B i * 
j =1 J J - 


1 T 


T J T 
,2 


T J 


l j 

'1 

4 


(13) 
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and thus the matrix P(k)H(k) required in the update of BOO, 
and equal to L(k) in Equation ( 1 1 ) , may be approximated by 


P(k)H(k) s x (k)[z(k) - y (k ) ] 



(14) 


where 

E[v?]= E[v 2 ] = ( 1 /2) a 2 ; E[z 2 ( j )+y 2 ( j ) ] = A +o y . 

in Equation (14) the vector X 0<) is <* ata independent and thus 
could be precomputed for k = 1,2,...,N for computational 

simplification of Equation (11). Similarly, in the implemen- 
tation of Equation (10), the first matrix may be replaced by 
the data independent matrix 


( .I^ N " j Bj)(A 2 +o 2 )T|. 

J = 1 J 


3(a). Modified Least Squares Algorithm 

If the noise £(k) in the signal model of Equation (8) 
were white, then the estimate B(N) obtained from the algorithm 
(10) or (11) would approach 0 as N - if one ignores the 
approximation made in arriving at Equation (8) and X is se- 
lected equal to 1. However, as the noise £(k) in Equation (8) 
is colored, there would be considerable bias in the parameter 
estimates under low to medium signal-to-noise ratios. To 
reduce such a bias or possibly eliminate it, the following 
simple modification is proposed. If the instantaneous fre- 
quency to ( t k ) at time t r given by (w 0 +2 W 36 ©*^ appearing in 
Equation (7) is small compared to 1/T S , then the noise vector 
£( k) is equal to v(k)-v(k-1) where v (k) i [v^k) v q (k)]'. 

This situation is illustrated in FIG. 2(a). 
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To eliminate the bias, the noise £(k) must be whitened 
by passing through the transfer function (l-z -1 )” 1 as shown in 

A 

FIG. 2(b) where denotes an unbiased estimate of 3 . The 

least squares algorithm is, in general, nonlinear and time- 
5 varying. However, if it is assumed that the algorithm in 

Equation (11) asymptotically approaches a time-i nvariant 
system, then under such an assumption, one may interchange the 
least squares algorithm with the transfer function (1-z”^)"~^ 
to arrive at the arrangement of FIG. 2(c). This, of course, 
10 corresponds to post-averaging the least squares estimates, 
which is the preferred technique for the present invention. 
Such a simple procedure provides very significant improvement 
in the estimates 3 ( k ) when the signal-to-noise ratio is low. 
In the simulations of the next section, the infinite time 
15 averaging (1-z -1 )*" 1 of FIG. 2(c) is replaced by an exponen- 
tially data-weighted averaging to take into account the time 
variation of the parameters to be estimated. 

3(b). DLS Algorithm in the Data Modulation Case 
20 In this case, the data samples D(k) in the signal model 

of Equation (4) take only possible values ±1 and the received 
signal may equivalently be written as: 

y (k ) = D(k)ASin( 0 (k))+v i (k) 

25 z(k) - D(k ) ACos (G( k ) ) +Vq(k ) ;k=1,2,...,N (15) 

Thus, as may be easily verified over any bit interval T b , 
where D(k) remains constant, the differential signal model of 
Equation (6) remains valid irrespective of the value of D(k). 
30 The model of Equation (6), however, is not applicable for 
those samples which lie on the bit boundaries, i.e., when y(k) 
and y(k-1) lie in different bit intervals. A simple modifica- 
tion of the algorithm to take care of the data modulation case 
is to simply discard such differential samples. If the number 
35 of samples M over any bit period is fairly large, this would 
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incur a negligible loss in the effective signal-to-noise ratio 
compared to the case of no data modulation. In fact, such a 
loss is simply equal to 1 01og 1Q ( 1-1/M)dB which is O.MB dB for 
M=1 0. This is corroborated by the simulations of the next 
section. 


3(c) Estimation of Time-Varying Parameters 

In the signal model considered above, it is assumed 
that the input signal parameter vector ip I0 is either a con- 
stant or a slowly varying function of time. In practice, this 
may be the case only over relatively short intervals of time, 
but there may be large variations in \Jijq over a comparatively 
large observation period. To take into account such a varia- 
tion and to ensure that the instantaneous difference frequency 
Q(t) f d/dt(0(t) ) (the sampled version of 0(t) given in Equa- 
tion (4)) remains within the low-pass filter-pass band of FIG. 
1, the parameter vector 4^0 generating the instantaneous fre- 
quency of the VCO is updated at regular intervals of T NT 3 sec 
for some integer N. The parameter vectors ijjp <|< Lo and i|» 0 
would change their values at intervals of T sec, assuming that 
the value of N is selected to be sufficiently small so that 
the variation in ^ over any T sec interval is small. Denot- 
ing by 0 Lq (T+) ,w Lo (T+) , etc., the values of reference oscilla- 
tor parameters just after the update at time T, are: 


“Lo (T+) = «Lo(T-)^ o (0/T).2T; o ( 0 / t )-3T 2 6 o (0/T) 

y Lo< T+) = r Lo (T-) + 3T6 Q (0/T) 

5 Lo (T+) ' A Lo (T_)+ So (0/T) (16) 

In Equation (16), 0 Lo (T-), fl Lo (T-), etc -’ represent the oscil- 

lator instantaneous phase, frequency, etc., just bei ore the 
correction, and the remaining terms on the right hand side 
represent the correction made on the basis of the estimation 
algorithm. Thus, 
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0 Lo (t-) = 9 Lo (0 + ) + , Lo (0 + )t + y Lo (0Ot 2 + « Lo (0 + )t 3 

fi Lo (T-) - u) L o (0 + ) + 2 Y Lo (0 + )T + 36 Lo (0 + )T 2 

r(T-) = y L o (° + ) + 3T6 Lo (° + ) 

= 5 lo< 0+) (17) 

Note that in Equation ( 1 6 ) there is no correction in the 
oscillator phase as the DLS algorithm does not provide the 
phase estimate. In Equation (16) <3 o (0+/T) denotes the esti- 
mate of parameter u> o (0+) obtained on the basis of measurements 
up to time T. Since there is no step change in the oscillator 
phase, the sampled measurements y(N) and z(N) at the demodula- 
tor output are the same with or without a correction at the 
instance NT S . However, the subsequent measurements y(N+j) and 
z(N+j) are now expressed with respect to the new parameter 
vector ij; 0 (T+)=[e o (T+) w 0 (T+) Y 0 (T+) 6 0 (T+)]* f ho (T+) ~^Lo (T+) 
as in Equation (18) below. 

y(N+j) = ASin(0(N+J))+v i (N+j) 

z(N+j) = ACos (0(N+J)) + Vq(N+j) ;j=1 N 

0(N+j ) = 0 o (T + )+co o (T+)jT s +T o (T+)(jT s ) 2 +6 o (T+)(jT s ) 2 (18) 

The last three elements of the vector i(j 0 (T+) will be zero if 
there is no change in the input signal parameters over the T 
sec interval and the estimate of ip o (0+) is obtained with zero 
estimation error. Thus, the technique is to set the a-priori 
estimate of the vector 4 i q (T+) equal to 0 and apply the DLS 
algorithm to estimate ip 0 ( T+ ) on the basis of observations 
(y (N+j ) , z ( N+ j ) ; j = 1 N } . The measurement model is ob- 

tained by simply replacing the index k by k+N in 
y(k),z(k),5 i (k),5 q (k) in Equations (6-9) but with i k =(k-1/2)T s 
as before (corresponding to a shift in time reference) . 

In the estimation of ^ 0 (T+) via the recursive algorithm 
of Equation (11) with the index k=N+1 , . . . , 2N , the "initial" 
covariance matrix P(N+1) is obtained as: 
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P( N+1 ) = AFP(N)F'+Q 


(19) 


where 


5 



*1 2T 3T 2 
0 1 3T 

-° 1 h 


and the matrix Q represents the uncertainty introduced due to 
the change in the input process parameters over the interval 
10 of T sec. Specifically, the last diagonal element of Q repre- 
sents the variance of the change in the parameter 66 (equal to 
the second derivative of frequency and related to the jerk of 
the physical trajectory) over the interval T. The above 
procedure is then extended in a straightf orward manner to the 
15 subsequent update intervals. The estimates of input signal 
phase and frequency at time instances £T+ are then simply 
given by e LQ UT+) and id Lo (!,T+) respectively for £=0,1,2 


4. Multistage Estimation 

20 Most of the phase and frequency estimation schemes for 

the first stage of a multistage estimate can be represented as 
in FIG. 3(a) in which elements common with FIG. 1 are identi- 
fied by the same reference numerals. An update algorithm 
(represented by a functional block 21 ) for generating the 
25 correction signal for the VCO could be any recursive or 
semi-recursive algorithm including an EKF or DLS algorithm, 
and the VCO update interval may be some integer multiple of 
the sampling period T s . FIG. 3(b) illustrates a functional 
block model for the generalized first-stage estimator illus- 
30 trated in FIG. 3(a). This functional block model will be used 
hereinafter as a symbol to illustrate the first stage in a 
multistage estimator shown in FIG. 4 in accordance with the 
objective of this invention. Note that 0(t), v.(t) and v (t) 

of FIG. 3(a) are respectively equal to G^t), v|(t), vht) of 



F89/31 4 


18 


FIG. 3(b). FIG. 3(c) illustrates a functional block model for 
subsequent stages of the multistage estimator of FIG. 4. 

An important observation made from FIG. 3(b) is that 
along with the estimate 0 j ( t ) of the input phase process 
5 Oj(t), there is also available a pair of signals, yVt) = 
AS i nO 4 1 ( t ) + vj(t), z^t) = ACosG(t) + v^(t) dependent upon 
the estimation error 0(t) = OjUJ-e^t). These error signals 
have exactly the same form as the signals at the input to the 
estimator. Moreover, the additive noise associated with these 
10 signals has statistics identical with the statistics of the 
noise at the input to the estimator. Therefore, this leads to 
the interesting possibility of the present invention, namely 
estimating the error signals 0(t) in a way similar to the 
estimation of ©.(t). In fact, the procedure can be repeated 
15 any number of times as shown in FIG. 4, although only a two- 
stage estimator is described below as an example of the pres- 
ent invention. It should also be noted that only the first 
estimation stage shown in FIG. 3(a) requires a VCO for down 
conversion because the actual input to this stage is at rf 
20 frequency u> c , whereas subsequent estimator stages generate the 
error signals by simple baseband computations. For example, in 
the di scr et e— ti me version of the estimation procedure, the 
signal at the output of estimation stage m may simply be 
computed as 
25 

y^k) = y m_1 (k)Cos(0 m-1 (k-1 ))-z m-1 (k)Sin(0 m-1 (k-1) ) 

( 20 ) 

z m (k)=y m ' 1 (k)Sin(0 m ' 1 (k-1)) + z m " 1 (k)Cos(0 ,n ' 1 (k-1)) ;m=2,3, ...n 

30 The refined estimate of the phase process 0j(t) in FIG. 4 at 
sampling instance k is then simply given by 


( 21 ) 


0 l(k) = 0 L (k)+ Gl (k)+.. . + 6n (k ) 
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The advantage of such a recursive estimation procedure is that 
the overall threshold in terms of the required CNR for the 
multistage estimator can be made much smaller than a single 
stage estimator, especially in situations involving high 
5 dynamics. 

In the prior-art single-stage estimator shown in FIG. 
1, the process parameters may be assumed to remain constant 
only over short intervals of time due to the high dynamics 
involved in a receiver carried, for example, in a high per- 
10 formance aircraft or spacecraft. Thus, the estimator is 
forced to use a relatively large noise bandwidth (shorter 
averaging period), resulting in large errors in the phase 
and/or frequency estimates. If the estimation errors are 
outside the region over which the error model (linear) assumed 
15 for the estimator remains valid, the estimator is said to be 
working below threshold or in the out-of-lock condition. In 
this condition, the estimation errors can be several orders of 
magnitude higher compared to the operation above threshold. 

In the multistage estimator of the present invention, 
20 this difficulty can be circumvented by successive reduction of 
the dynamics (the estimation errors due to dynamics) at the 
output of consecutive estimator stages and by averaging the 
signal over progressively longer intervals (and thus progres- 
sively reducing the effect of noise) over which the process 
25 parameters remain nearly constant. In this multistage estima- 
tion structure, none of the individual stages (except the last 
one) need necessarily operate above its threshold. For the 
convergence of the overall estimator, one only requires that 
the estimates are made in the right direction (estimation 
30 errors do not exceed the parameters to be estimated in some 
average sense). 

5. Recursive DLS-EKF Algorithm 

The following section considers a simple special case 
of n=2 for the more general multistage estimator of FIG. M, 


35 
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wherein the first stage is as described above with reference 
to FIGs. 3a and 3b and the DLS algorithm is as described in 
the previous section above. The second stage is an extended 
Kalman filter (EKF ) as described above with reference to FIG. 

5 3c. As the dynamic variation of the error signal 0(k) - 0 (k) 

of the first stage is much smaller compared to the original 
signal 0(k) (the frequency variation over any update interval 
is much smaller), the effective averaging time period for the 
second-stage estimator implemented with a Kalman filter can be 
10 selected to be higher than for the first stage implemented 
with a DLS algorithm. This is achieved by selecting a smaller 
value of the "dynamic noise" covariance matrix Q for the EKF 
and/or a higher value for the exponential data weighting 
coefficient X of Equation (9). 

15 

5a. Extended Kalman Filter (EKF) Algorithm 

Referring to FIG. M, consider the problem of estimating 

the unknown error signal parameters anc * ^oZ 

Zth VCO update period for any integer 121 on the basis of the 
20 set of measurements {y 1 (k) .z 1 (k) } of Equation (22) below by an 

EKF 


25 


30 


35 


y 1 (k) 
z?(k) 
e^k) 


= ASin(0 1 (k))+v- (k) 


ACos(0 1 (k))+Vq(k) 
k = N( Z-1 ) + j ; j=1 » 2, • • • > N ;Z=1,2, 


( 22 ) 


Note that as for the first-stage DLS algorithm the parameter 

, r r\ „ Y * „ V may be different over differ- 

vector i|; 0 £ = [00^ ^ 0 Z 'ol 6 oZ J y 

ent VCO update intervals. For computational simplicity, a 
third-order EKF is used in the second stage and the contribu- 
tion of the last term In the expression for £>’ 00 Is Ignored, 
which is appropriate for the GPS trajectories considered here. 

Denoting the state and parameter vectors at time k = NU-D+J 
by ^(j) and n* respectively, i.e., with «|»(j) = H J T s 



F89/31 4 


21 


0.5(jT s ) 2 ]', n 4 = [e 0 £ w 0 £ 2Y 0 £]’, the extended Kalman 

filter equations for the update of n 4 in the second stage, the 
estimate of n 4 » are given by 

5 = fi 4 (j-1 )+K 4 (j )v £ (j ) 

K 4 ( j) = Z A CJ“i)*CJ >) -1 

I 4 (j) = ll£(j-1)"I ]l (j-l)<P(j)[X+^ , (j)I 4 (j-1)<|<(j)]" 1 

♦ •(JJlitj-DJA (23) 

10 v 4 (j) = y 1 (k)Cos(<j> 4 (j))-z 1 (k)Sin((j) 4 (j)) 

4> 4 (j) = ip' (j )h 4 ( J-1 ) ; k=( A-1 )N+j; j=1 , 2, . . . ,N;A=1 , 2, . . . 

In the Equations (23) above, the initial estimate fi 4 (0) is 
simply taken to be equal to n 4 _-j(N). This is an appropriate 
15 choice for the initial estimate in view of the fact that if 
the first stage of the DLS estimation algorithm is convergent 
then n 4 would possess some continuous drift term, i.e., w o4 
will have a component linear in time if Y o4 ^ 0. The "initial 
error covariance" matrix £ 4 (0) is simply set equal to some 
20 diagonal matrix representing the uncertainty in the difference 
parameter n£ - n 4 _i» 

5b. Estimation in the Presence of Data Modulation 

In this case one could apply the more sophisticated 
25 version of R. Kumar (1988) cited above, wherein an explicit 
detection of possible data transitions is followed by the 
demodulation of data, thus effectively reducing the problem to 
the case of no data modulation considered above. However, 
here such a detection is bypassed, and instead a simple modi- 
30 fication in the estimation algorithm is proposed that takes 
into account the data modulation. If the VCO update interval 
T is selected equal to bit period T b , then the data modulation 
represents an additional phase uncertainty at the boundaries 
of the update intervals. This is taken into account by adding 
an appropriate value, say (tt/2) 2 , to the first diagonal ele- 
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10 


15 


merit of ^(0) and modifying the initial estimate fj^(0) by tt/2, 
i.e., n^(0) = (N) +tt/2, for those values of i that corre- 

spond to bit boundaries* Such an algorithm is expected to 

result in somewhat higher estimation errors compared to the 
more sophisticated scheme of R. Kumar 0 988) cited above, but 
is much simpler in terms of implementation. In the case of a 
two-stage estimator (stages 1 and 2 of FIG. 4) , the estimates 
of the input signal phase and frequency at time instances IT 
are given by 


OUT) - 0 L (U- 1 )T+) (N)fi 5L (N) 

OUT) = 0 l (U~1)T+) + <3 0 £ + 2Y Q £ * ( NT S ) 

where and represent the second and third 

respectively of rj^(N). 


(24) 


elem ent 


6. Simulations 

Simulation results obtained when the algorithm is 
applied to the tracking of phase and frequency for high dy- 
namic GPS receivers will now be presented. For the purposes 
20 of simulation, assume that the pseudo-random code has been 
removed from the received signal, and symbol timing has been 
acquired. The simulation considers both the case of the data 
modulation removed via an auxiliary link and the case of an 
unknown data modulation present. For the simulations a sam- 
25 pling rate of 500 samples/second is assumed with a high dy- 
namic trajectory considered previously in V.A. Vilnrotter, et 
al . , "A Comparison of Frequency Estimation Techniques for High 
Dynamic Trajectories," JPL 88-21, Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena, CA, September 
30 1988, and reproduced in FIG. 5 where the acceleration and the 

jerk (the derivative of acceleration) are measured in units of 
g (the gravitational constant equal to 9.8m/s). In the case 
when data modulation is present, a BPSK modulation at a rate 
of 50 bits/second is assumed. 
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The parameters of most interest in this application are 
the instantaneous phase and frequency of the input signal 
r^Ct), which corresponds to the high dynamic GPS trajectory of 
FIG. 5. Since the present invention is mainly interested in 
the tracking performance of the proposed algorithm, it is 
assumed that the initial trajectory parameters at zero time 
are known as in V. Vilnrotter, et al., Sept. 1988, supra. The 
received signal carrier frequency f Q = u c /2tt in the signal 
model of Equation (1) is taken to be equal to 1.575 GHz. The 
GPS receiver instantaneous pseudo-range R in meters and veloc- 
ity v d in m/s are related to the instantaneous phase 0-[(t) of 
Equation (1) and its derivative (5 j ( t ) as 


Q c 


v d = f d f c = 2 tt f c 

where f d denotes the Doppler in Hz and c is the speed of 
light. Denoting by RU) and ? d U) the estimates for RU) and 
f d ( ^ ) res P ec ki vely , which denote the range and Doppler of the 
input trajectory at the end of the &th update interval, then 
the performance measures of the estimation algorithm are given 
by the following sample rms values of the estimation errors 



L 

£[rU)-rU)] 2 

£=i 


d,rms * 


A 


£[f d U)-f d U)y 
8. =1 


( 26 ) 


where L = 4000/N is the number of update intervals for the 
entire trajectory. These measures are obtained as a function 
of P/N q , where P denotes the received carrier power and N Q is 
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the one-sided power spectral density of the receiver bandpass 
noise. 

At lower range of (P/N Q ) ratio, the receiver may lose 
frequency lock in that the frequency errors at times may 
5 exceed ± one-half of the sampling frequency, f or ± 250 Hz. 

Since the error signals e(k) of Equation (11) are the same for 
frequency errors of A Hz as for the case of A + nf s Hz for any 
signed integer n, the estimator may make frequency estimation 
errors of nf s Hz. This situation may be referred to as cycle 
10 slipping in the frequency estimator and is akin to the phe- 

nomenon of cycle slipping (phase errors equal to multiples of 
2 tt) in the phase-estimators. If there are one or more cycles 
slipped in frequency, the computed value of f d pm3 would be 
much larger compared to the case when no such cycle slips 
15 occur and would be unacceptable. Thus, another important 

parameter for the performance is the probability of maintain- 
ing frequency lock throughout the trajectory denoted P(lock) 
or the probability of losing the lock P^=1 -P(lock) . For the 
purposes of estimating the probability by digital computer 
20 simulations, 100 simulation runs are made for each value of 

P/N q of interest and an estimate of P^ is plotted vs the 
carrier power to noise power spectral density ratio (CNR). The 
sample rms values of Equation (26) are also averaged over all 
those simulation runs for which the frequency lock is main- 
25 tained. It may well be that for sequences that remain under 
frequency lock there may be slipping of cycles in the phase 
estimates. However, even under the presence of such cycle 
slips, the computation made on the basis of Equation (25) 
provides a good estimate of the pseudo-range as evidenced by 
30 the simulations. One cycle slip only corresponds to an error 
of about 0.2 meters in the pseudo-range estimate. 

FIGs. 6-15 present the simulation results for' the DLS 
algorithm and the composite DLS-EKF algorithms presented in 
the previous sections. The results for the EKF algorithm 
35 operating by itself are available in V. Vilnrotter, et al . , 
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Sept. 1988, supra, for comparison. For the simulation results 
a value of A equal to 0.97 has been selected. The initial 
covariance matrix P(0) for the DLS algorithm is selected to be 
a diagonal matrix with its diagonal elements equal to 2x10^, 
2x10' and 2x10^ respectively, reflecting the possible uncer- 
tainty about the parameters. Three different VCO update 
intervals equal to 5, 10 and 20 sample times have been consid- 
ered. The matrix Q of Equation (19) is also selected to be a 
diagonal matrix for convenience, with its consecutive diagonal 
elements equal to 4x10^, 2x10® and 10®. The Q matrix repre- 
sents possible variations in the input signal parameters over 
an update interval and is arrived at from the consideration of 
a-priori estimate of the maximum possible value of the highest 
order derivative (jerk) present in the input trajectory. For 
the second stage EKF algorithm, the initial covariance matrix 
X 4 (0) is selected to be also a diagonal matrix but with its 
elements smaller in value than the corresponding elements of 
the Q matrix, thus effectively resulting in a higher averaging 
period and smaller estimation errors compared to the DLS 
algorithm. The selected values of diagonal elements of 1^(0) 
matrix are equal to 1.0, 10^ and 10® respectively, in the 
following simulations. From the simulations it appeared to be 
advantageous in terms of numerical stability to periodically 
reset the covariance matrix P of the DLS algorithm to its 
initial value. Such a period was selected to be 10 times the 
VCO update interval. 

FIGs. 6 and 7 present the simulation results for the 
performance of DLS algorithm while tracking the high dynamic 
trajectory of FIG. 5 in the absence of any data modulation. 
FIG. 6 plots the probability of losing the frequency lock 
as a function of CNR for two different values of N equal to 10 
and 20. As may be inferred from the figure, a value of p^ of 
less than 0.1 is obtained for CNRs above 23.1 dB which is 
defined to be the threshold point of the algorithm. FIG. 6 
also plots the average number of cycles slipped in the fre- 
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quency estimation, denoted by N cs . In defining such a cycle 
slip, the entire frequency range is divided into disjoint 
segments of f g Hz with the first segment extending from -f g /2 
to f /2 Hz. Whenever the frequency estimation error jumps 

o 

5 from one such segment to an adjacent one in either direction, 
a cycle slip is said to occur. FIG. 7 plots the rms error in 
the Doppler estimation as computed from Equation (26) and 
averaged over all convergent sequences. For a CNR between 25 
and 30 dB-Hz , an rms error of 10-20 Hz is obtained that corre- 
10 sponds to a velocity tracking error of 2-4 m/s. FIG. 7 also 
plots the average length L cs of a slipped cycle in terms of 
number of samples. The information about N cs and L cs is 
relevant in the case of multistage algorithm. FIGs. 8 and 9 
present the results for the probability of losing lock and 
15 the rms estimation error for the DLS algorithm in the presence 
of data modulation for three different values of N equal to 5, 
10 and 20. As may be observed from the figures, the presence 
of data modulation increases the threshold by only 0.25*0.5 dB 
compared to the case of no data modulation. The increase in 
20 rms frequency estimation error is about 1 0% due to data modu- 
lation. 

FIGs. 10 and 11 present the performance of the compos- 
ite DLS-EKF algorithm in the absence of data modulation. Note 
that corresponding to N=5, the threshold of the algorithm is 
25 22.75 dB-Hz, which is slightly smaller than for the DLS algo- 

rithm. However, the rms estimation errors are significantly 
smaller than for the single-stage DLS algorithm. For the CNR 
range of 25*30 dB-Hz, the rms error in the Doppler estimation 
lies in the range of 4-15 Hz corresponding to the velocity 
30 estimation error range of 0.8 to 3 m/s. The DLS-EKF algorithm 
also provides the carrier phase estimate. The modulo-2n 
phase- estimation error is plotted in FIG. 12 from which it is 
clear that the algorithm is capable of coherent data detection 
with small probability of error if the CNR is higher than 25 
35 dB-Hz. In fact, as shown in FIG. 15, the algorithm provides 
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good estimates of pseudo-range (related to the absolute phase 
error via Equation (25)) up to a CNR of about 2 3 dB at which 
point the rms error is less than M m. The rms pseudo range 
error is less than 1 m for CNRs higher than 25*5 dB-Hz. 

5 The corr es ponding results for the performance of 

DLS-EKF algorithm in the presence of data modulation are 
presented in FIGs. 13~15. For this case, a minimum threshold 
of 23.8 dB is obtained for N=10 which is about 1 dB higher 
than for the case of no data modulation. In terms of rms 
10 estimation errors , for a CNR range of about 25-30 dB-Hz, the 
rms frequency estimation error lies in a range of 8-20 Hz 
corr esponding to a velocity error of about 1.5 to ^ m/s. For 
this case, as is apparent from FIG. 15, the pseudo-range 

estimation errors are also higher and for a CNR range of 25-30 
15 dB-Hz lie in a range of 0.3-6 m. Notice, however, that no 

sharp threshold is observed in either the frequency or phase 

estimation errors over the entire range of CNR between 22-30 
dB-Hz considered in the simulations. 

20 7. Comparison with Previous Techniques 

For the case of no data modulation, the following 

compares the performance of the proposed algorithm with some 
of the techniques analyzed in V. Vilnrotter, et al . , (1 988) 

cited above, in terms of their performance when tracking 
25 exactly the same high dynamic trajectory. Compared to a more 
computation-intensive maximum likelihood estimate, the DLS-EKF 
algorithm requires about 0.25 dB smaller CNR than MLE in terms 
of threshold. In terms of rms frequency estimation errors, 
the MLE achieves an rms error between 8 Hz to 35 Hz at a CNR 
30 of 23 dB-Hz depending upon the estimation delay ranging be- 
tween 30-80 samples (higher delay provides smaller error). The 
DLS-EKF algorithm provides an error of 35 Hz for a delay of 5 
samples at a CNR of 23 dB-Hz. The MLE algorithm does not 
provide any phase estimate. 
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25 


compared to a single-stage EKF algorithm, the DLS-EKF 
algorithm is about 1. 5-2.0 dB better in terms of threshold 
depending upon the value of exponential data weighting coeffi- 
cient and the filter order used in the simulations of V. 
Vilnrotter, et al . , (1988) cited above. In terms of rms er- 
rors, the performance is similar to that of third-order EKF 
alone. Notice, however, that direct comparison with the 
results of V. Vilnrotter, et al . , may be somewhat misleading. 
This is so because while the DLS-EKF algorithm includes all of 
the sequences in the computation of rms error above 25.5 
dB-Hz , EKF rejects about 5* of the worst sequences, as the 
probability of losing loch is about .05 at CNR of 25.5 dB-Hz. 

For the case of cross product AFC loop analyzed by V. 
Vilnrotter, et al., the threshold lies in a range of 25-28 
dB-Hz depending upon the loop parameters. Thus, the DLS EKF 
algorithm is superior by 2-5 dB-Hz compared to AFC loop. AFC 
loop provides a minimum rms frequency error of 25 Hz at a CNR 
of 28 dB-Hz compared to a minimum of 5 Hz achieved 
DLS-EKF algorithm for the same CNR. Notice that in an AFC 
loop, the parameters achieving a relatively low estimation 
error are different than those yielding low thresholds and 
thus, a range of loop parameters must be considered for proper 
comparison. In terms of rms phase error, the performance of 
the DLS-EKF algorithm is similar to EKF alone. In 
computations, the DLS-EKF algorithm requires about the same 
number of computations as for a fourth-order EKF but about 
twice as many computations as a third-order EKF. The number of 
computations are at least an order of magnitude smaller than 

the MLE. 

For the case when the data modulation is present, the 
following compares the performance of the DLS-EKF algorithm 
with the MLE algorithm of W.J. Hurd, et al . , (1980 

above, where a somewhat less severe CPS trajectory is ana- 
lyzed. The results of W.J. Hurd, et al . , show a marked 
35 threshold of about 30 dB-Hz in terms of CNR compared to less 
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than 24 dB-Hz for the proposed algorithm. Thus, the proposed 
algorithm results in more than 6-dB reduction in terms of 
threshold compared to previous schemes of the literature. 

In terms of the rms frequency estimation errors at a 30 
dB-Hz CNR, the MLE algorithm provides an rms range rate error 
of 6 m/s compared to an error of less than 2 m/s achieved by 
DLS-EKF algorithm. There is also very significant improvement 
in terms of the rms position estimation error. At 30 dB-Hz an 
rms error of 1 meter is reported by W.S. Hurd, et al . , com- 
pared to about 0.25 meter obtained by the DLS-EKF algorithm. 
In terms of computations, both of the algorithms are compara- 
ble. In terms of threshold on CNR, the DLS algorithm is very 
close to the composite DLS-EKF algorithm. However, in terms 
of rms frequency estimation errors , it has significantly 
higher estimation errors.' In those cases where higher estima- 
tion errors are acceptable, one may apply the DLS algorithm by 
itself, as it requires only one-half of the computations 
required by the DLS-EKF algorithm. 

8. Conclusions 

A novel multistage estimation scheme has been presented 
for the efficient estimation of the phase and frequency of a 
very high dynamic signal, which may possibly be phase modu- 
lated by unknown binary data and is received under relatively 
low carri er-to-noise power ratio conditions. The proposed 
scheme is of very general nature and has much wider scope than 
the applications described as examples. For a very important 
application of dynamic GPS trajectories, the specific example 
of a DLS-EKF scheme described has been specialized to have 
just two stages. The first stage of the estimation scheme is 
a least-squares algorithm operating upon the differential 
signal model while the second 3 tage is an extended Kalman 
filter of third order. 

For very high dynamic GPS traj ectori es , the proposed 
algorithm has been shown to significantly outperform the 
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previous algorithms reported in the literature in one or more 
aspects, including threshold on CNR, estimation errors, avail- 
ability of phase estimates and thus the estimate of pseudo- 
range, computational complexity and flexibility. For the case 
of no data modulation, the proposed scheme has a threshold 
that is slightly lower than the more computation-intensive 
implementation of MLE algorithm. When compared to just the 
EKF operating by itself, the proposed DLS-EKF scheme provides 
from about 1.5 to 2-dB reduction in threshold. In comparison 
to more conventional schemes, such as AFC loops, the perform- 
ance is even better. 

For the case when an unknown data modulation is pres- 
ent, the algorithm provides an improvement of 6 dB in terms of 
threshold on CNR in comparison to the MLE scheme of W.J. Hurd, 
et al . , (1 987) specifically proposed for such applications. 

In addition to phase and frequency estimates, the algorithm 
can provide estimates of frequency derivative as well, al- 
though not presented here. The scheme being of a very general 
nature, it may be possible to reduce the threshold even fur- 
ther by using a higher dimension for the state vector related 
to the higher number of terms in the Taylor series expansion 
in arriving at the signal model for the first stage DL^> algo 
rithm. Further improvements are possible by the application 
of more optimum sampling techniques as proposed by R. Kumar, 
in "Differential Sampling for Fast Frequency Acquisition via 
Adaptive Extended Least Squares Algorithm," (1987) cited 
above. The performance may also be improved both in terms of 
the threshold and the rms estimation errors by increasing the 
number of states to three or more. 
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MULTISTAGE ESTIMATION OF RECEIVED 
CARRIER SIGNAL PARAMETERS UNDER VERT 
HIGH DYNAMIC CONDITIONS OF THE RECEIVER 


ABSTRACT OF THE DISCLOSURE 

A multistage estimator is provided for the parameters 
of a received carrier signal possibly phase-moduLated by 
5 unknown data and experiencing very high Doppler, Doppler rate, 
etc., as may arise, for example, in the case of Global Posi- 
tioning Systems (GPS) where the signal parameters are directly 
related to the position, velocity and jerk of the GPS ground- 
based receiver. In a two-stage embodiment of the more general 
10 multistage scheme, the first stage, selected to be a modified 
least squares algorithm referred to as differential least 
squares (DLS), operates as a coarse estimator resulting in 
higher rms estimation errors but with a relatively small 
probability of the frequency estimation error exceeding one- 
15 half of the sampling frequency, provides relatively coarse 
estimates of the frequency and its derivatives. The second 
stage of the estimator, an extended Kalman filter (EKF), 
operates on the error signal available from the first stage 
refining the overall estimates of the phase along with a more 
20 refined estimate of frequency as well and in the process also 
reduces the number of cycle slips. 
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